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Abstract 

Quantum phase transitions and observables of interest of the ground state in the Tavis- 
Cummings model are analyzed, for any number of atoms, by using a tensorial product of 
coherent states. It is found that this "trial" state constitutes a very good approximation to 
the exact quantum solution, in that it globally reproduces the expectation values of the matter 
and field observables. These include the population and dipole moments of the two-level atoms 
and the squeezing parameter. Agreement in the field-matter entanglement and in the fidelity 
measures, of interest in quantum information theory, is also found. The analysis is carried out 
in all three regions defined by the separatrix which gives rise to the quantum phase transitions. 
It is argued that this agreement is due to the gaussian structure of the probability distributions 
of the constant of motion and the number of photons. The expectation values of the ground 
state observables are given in analytic form, and the change of the ground state structure of 
the system when the separatrix is crossed is also studied. 



PACS numbers: 42.50.Ct, 64.70.Tg, 03.65.Fd 



I. INTRODUCTION 



A model of fundamental importance in quantum optics and essential in any descrip- 
tion of systems involving the^interaction between matter and light is given by the Hamil- 
tonian proposed by Dicke [1!]. This describes the interaction of N-identical two- level 
systems with a single mode radiation field in the dipole approximation; the particles are 
confined in a container small compared to the radiation wavelength. 

The simplest completely soluble quantum-mechanical model of one atom in an elec- 
tromagnetic field is described by the Jaynes-Cummings (JC) model In this 
model the Hamiltonian terms which do not conserve the energy of the field plus the 
occupation number of the two-level system have been neglected because their contribu- 
tions are very small except for very high intensity fields which has been called the 
rotating-wave approximation. 

Since the end of the sixties, the exact solution for an N-molecule - radiation field 
Hamiltonian was obtained [5|, and a full family of exactly solvable generalizations 
has been presented recently |6|. They have provided insight into other more compli- 
cated physical systems and offers a standard comparison for approximation techniques. 
Presently, the JC model and its generalization given by Tavis and Cummings, the Tavis- 
Cummings (TC) model, continue to be fundamental to study basic properties of quantum 
electrodynamics and to understand phenomena like the existence of collapse and revivals 
in the Rabi oscillations \l\ , the formation of macroscopic quantum states, and the many 
experimental studies of Ry^dberg atoms with very large principal quantum number within 
single-mode cavities |3|, U, Isj. It is important to mention that it was not until 1987 that 
the collapse and revival predicted by the JC model was experimentally observed j^. 

Instabilities, or quantum phase transitions, in these models have recently been stud- 
ied due to their influence on several kinds of entanglement, important in quantum in- 



formation theory. Lambert et al. 



10| have studied the atom-field entanglement in the 



thermodynamic limit, finding logarithmic divergences in the atom-field entanglement 
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and discontinuities in the average linear entropy. Buzek et al. [llj showed that an in- 
finite sequence of instabihties appear in the ground state of the Dicke model when the 
rotating-wave approximation is taken, characterized by changes in the bipartite entan- 
glement between atoms. Similar models have been used to describe the phase transitions 



of identical qubits interacting with a bosonic mode in an adiabatic approximation 121]. 



Although there have been recent contributions to the quantum phase transitions in 



the Dicke model and its extensions 



3, 



14 



15|, we want to enhance the role of the 



catastrophe formalism to determine when significant changes to the ground state occur 
for small changes of the external environment (the parameters of the model). According 



to Sachdev [l6[, if a system undergoes a ground state energy phase transition at T = 
as a function of external parameters, then it also has a thermodynamic phase transition 
for fixed interaction parameters as a function of the increasing temperature. 

Quantum phase transitions and stability properties of algebraic models have been 
studied through the catastrophe formalism and the coherent states theory Recently, 
a procedure was established to determine the phase transitions associated to nuclear 
and matter models [is!]. In this contribution we present a comprehensive study of the 
phase transitions of the ground state in the TC model, for any number of atoms, when 
the strength of the interaction 7 and the atom hyperfine two-level separation energy 
uja are varied. These quantum transitions occur also for a small number of particles, 
although we can take the thermodynamic limit and corroborate the results obtained at 
the beginning of the 70's [l^. The minima of the expectation value of the Hamiltonian 
are calculated for a trial state, built as a coherent state HW{1) ® SU{2). Employing 
the catastrophe formalism, the separatrix is found. The structure of the ground state of 
the system changes when this separatrix, given by uja = ±7^, is crossed, in agreement 
with the phase transition found by Hepp and Lieb [19]. These parabolae divide the 
control parameter space in three regions which map naturally onto the Bloch sphere, 
where minima of the expectation value of the Hamiltonian with respect to the coherent 
trial state are attained: (1) the North Pole, where the ground state of the system is 



constituted by all the atoms lying in the low hyperfine level and the electromagnetic 
field has zero photons; (2) the South Pole, with the ground state formed by all the 
atoms in the excited hyperfine level and zero photons in the electromagnetic field; and 
(3) the Parallel region, where the ground state has a variable number of components 
(field + matter), ranging from one, for values of the interaction strength close to the 
arms of the parabolae, to N + 1, for the collective regime. By means of Ehrenfest's 
thermodynamic classification of phase transitions we have determined that for the TC 
model there are second order phase transitions when the separatrix is crossed through 
the arms of the parabolae, and first order transitions when the crossing happens through 
their vertices. We also explore the influence of the phase transitions on the behaviour of 
observables of interest for the matter and the field: the population and dipole moments 
of the two level atoms, together with the expectation value of the number of photons and 
its fiuctuations. We evaluate the squeezing and entanglement properties of the ground 
state in the different regions of the control parameter space. In order to compare the 
reduced probability distribution of the excited number of atoms, obtained by the exact 
quantum result, with that determined using the semi-classical approximation, the fidelity 
between them is calculated. Finally, the photon number distribution of the ground state 
as a function of the atom-field coupling parameter 7 is given. 



II. SEMICLASSICAL ANALYSIS 



By semiclassical analysis it is meant here the following procedure: i) calculation of the 
expectation value of the Hamiltonian with respect to the tensorial product of coherent 
states {trial state); ii) determination of the minimal critical points; and iii) use of the 
catastrophe formalism to find the stability properties. The use of coherent states as 
trial states lets us determine in analytical form the expectation values of matter and 
field observables. 
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A. Energy Surface and Critical Points 

We are interested in solving the TC Hamiltonian given by 



Htc = CL^CL + uj'j^Jz -\ — 1= (a^ + a J+) , (1) 

V N 



where is the number of particles in the system, uj'p is the field frequency, the atomic 
energy-level difference, and 7' the field-atom interaction strength. It is immediate that 
this Hamiltonian commutes with the operator 

A = a^a + J, . (2) 

It is then convenient to rewrite it by introducing a detunning parameter A, and by 
dividing it by cu^ (which can be thought of as the natural unit of frequency) and by the 
total number of particles, having in this way an intensive Hamiltonian operator 

where A = l = \ — uja and 7 = — . 

In order to obtain an energy surface we use as a trial state the direct product of 



coherent states in each subspace: Heisenberg-Weyl for the photon part [s] and SU{2) 



or spin for the particle part |20|, i.e., \a, () = \a) ® \0. Explicitly, it is given by the 
expression 

where the ket |z/) is an eigenstate of the photon number operator, |j, m) is a Dicke state 
with j denoting the eigenvalue of J^, and m the corresponding eigenvalue of Jz- The 
trial state contains particles and up to an infinite number of photons distributed into 
all the possible ways between the two levels. 
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The expectation value of the TC Hamiltonian ([sj]) in this trial state is defined as the 
energy surface given by 

E{q,p,e, 0) = {a,C\H\a,0 

+ A— cost^H sm tl [q cos (p — p sm cp) . [b) 



N N V2iV 

To get this expression we have substituted 

a = -^(g + ip) , (6) 
C = tan exp {i 0) , (7) 

where (p, q) correspond to the expectation values of the quadratures of the field and 
{9^ 0) determine a point on the Bloch sphere. The function A is the expectation value of 
the constant of motion A, and is given by 

A(j, g, p, 6, 0) = (a, C|A|a, C) = ^ + p') - 3 cos 9 . (8) 
The eigenvalue j is related to the number of molecules or atoms present in the two-level 

n 

system, and it can take the values \m\ < j < N/2 [20], according to the symmetry of the 
particles. If we are interested in the symmetric configuration, the so called superradiance 
regime |l9|, we have to consider N = 2j, which will be assumed hereafter. 

The critical points of the energy surface are obtained by equating its first derivatives 
to zero. The derivatives with respect to q and p give us the relation of them with the 
critical values of 9 and 0: 



Qc = -VJ7 sin 6'c COS0C , (9) 
Pc = a/7 7 sin 6'c sin 0c . (10) 



The derivative with respect to 9, using (9) and (10), determines its critical value 9c in 
terms of the strength parameter 7 and of cj^: 

sin 9c {uJ A — l"^ cos 9c) = . (11) 



The remaining equation, related with the derivative with respect to 0, gives us an 
identity. One can check that the Hessian matrix has rank 3 at the critical points. This 
means that the energy surfaces at the critical points are independent of (j) and we shall 
say the energy surface is 0-unstable. 

We thus obtain the following critical points: the North Pole {9c = 0), the South Pole 
{6c = vr), and a Parallel defined by 6c = arccos {ua/'y'^). 

We determine the nature of the critical points examining the hessian matrix of the 
energy surface at the critical points. Because the variable is spurious, we have the 
liberty to choose its value (which in general we will take to be zero). We can recover 
the values of quantities dependent on 0, at any (p, by applying the transformation e*"^^. 

In the North Pole {6c = 0) the hessian matrix has the three eigenvalues 

1/(2 j), (l + J^A ± V(1-j^a)2 + 4j7') /(4j) . 

For uja > 7^ we have a minimum (all three eigenvalues positive); for uja = 7^ we obtain a 
degenerate critical point; and for ua < 7^ we get a saddle point (one eigenvalue is always 
negative). Along tUa = 7^, the hessian matrix is singular and thus it is a bifurcation set 



of the system 17|. 



In the South Pole {6c = vr) the hessian matrix has eigenvalues 

1/(2 j), {±{1-Juja) + ^(r+T^T^F + ITy) /(4j) . 

Here, for ua < —7^ we have a minimum; for uja = —7^ we obtain a degenerate critical 
point; and for lja > —1^ we get a saddle point. The matrix is singular along oja = —7^, 
which is then another bifurcation set. 

Finally, in the Parallel region {6c = arccos(c(j^/7^)), the hessian matrix eigenvalues 

are 



1/(2 j), [l'{l + J 7') ± ^7^ (1 - J + 4 J 7^ j /(4 J 7') . 

In this case, there are no solutions above uja = 7^ nor below oja = —7^, and in between 
all the eigenvalues are positive. Along Ua = ±7^, the hessian matrix is singular. 
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The minimum energy function therefore maps all the region uoa > 7^ to the North 
Pole of the Bloch sphere, all the region uja < —7^ to the South Pole, and the region 
—7^ < < 7^ to the rest of the sphere. Along the curves ua = ±7^ the nature of 
the critical points is not determined by the hessian; they are bifurcation sets in the 
parameter space forming part of the separatrix of the system, which is shown in Fig.([l]). 



FIG. 1: On the left, the separatrix of the system in the parameter space (7,0;^) is displayed. 
A trajectory in the control parameter space is shown, along which we analyze the behaviour 
of the semiclassical energy surface. On the right, 9c at minima is shown as a function of the 
interaction strength 7 and uja- 

Table [T] shows the critical points, the energy, and the constant of motion evaluated 
at these, together with the conditions in the parameter space to guarantee that they 
constitute an energy minimum. 

Most system variables will inherit their behaviour from that of 6c at the minima. 
Fig.([l]) shows 6c as a function of the interaction strength 7 and uja- The full energy 
surface is shown in Fig. ([2]) as a function of 7 and uja also. 
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TABLE I: Expectation values of energy per number of particles and constant of motion, of 
the estimated ground state of the system. The last column shows the conditions for E to be a 
minimum in each region. 



Minima 


N 


A 


Conditions 


Oc = 


LOA 

2 


-j 


UJA > 7^ 




UA 
2 


j 


< -7^ 


6c = arccos (0^0/7^) 


472 


■ -u)A (<.^A+2)+7* 
J 272 


-72 < < 7^ 



B. Phase Transitions 



The order of the phase transitions can be determined following the Ehrenfest classi- 
fication. A phase transition takes place between p and q branches of critical points, and 
is of nth-order if 



lim 

5-^0 



(p), 
I 



for i 



0, 1, 2, 



sa-S 

n — 1, but fails for i 



lim 

5^0 



d'Ei;'\s) 



n 



(12) 



so+5 



17l |. In the TC model, the phase 
transitions for the states of mimimum energy occur at the separatrix uja = ±7^, as we 
have seen. It can be shown that crossing the arms of the parabolae, such as along the 
circuit I — II — III — IV — / in Fig.([T]), leads to second order phase transitions; in this 
case the composition of the ground state changes from only one component |0) ® |j, —j) 
at the North Pole to many components at the Parallels region, or from one component 
|0) ® \ji j) at the South Pole to many at the Parallels, and viceversa. Figure ([3]) below 
shows contour maps for E = E{6, (p) just before and after the crossing at 7 = — 1 along 
Path II. 

The crossing at the vertex 7 = u; = along a non-zero slope, such as Path V in 
Fig.([l]), is a first order phase transition. The minimum energy state switches from all 
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FIG. 2: Energy surface ii^ as a function of the interaction strength 7 and the atomic energy 
level difference loa- The plane growing with uja represents E at the North Pole, the other 
plane is E at the South Pole, and the paraboloid is E at the Parallels region. Note that the 
planes cut the paraboloid at \oja\ = 7^; thus, the minimum energy is found at different regions 
depending on the relationship of oja to 7. 

the atoms in their ground state and zero photons in the field (North Pole), to all the 
atoms in their hyperfine excited state and zero photons in the field (South Pole). This 
transition involves a change in sign for uja-, which may be achieved by continuously 
varying a magnetic field to which the atoms in the cavity are exposed. The energy 
surface E{9, (p) changes abruptly through this vertex, as expected, and it is shown in 
Fig.dl. 
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FIG. 3: Contour levels of the surface energy function E{9, (p) near the arms of the separatrix. 
Shown are the plots for 7 = —1.35 (left, Parallels region) and for 7 = —0.9 (right, North Pole), 
both at UJA = 1 and j = 1. The shape on the right remains fixed throughout the whole North 
Pole, and is distorted again only after the separatrix is crossed. Darker regions represent lower 
levels. 




6. 6- 



FIG. 4: Contour levels of the surface energy function E{9, </>) near the vertex of the separatrix. 
Shown are plots for 7 = loa = —0.001 (left. South Pole) and for 7 = loa = 0.001 (right, North 
Pole). The change is abrupt. Darker regions represent lower energy levels. 
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TABLE II: Expectation values of the quadratures of the field, photon number operator and 
their corresponding fluctuations, at the minima. 



Minima 


{q)/VN 




{h)/N 


(An)7iV2 


0c = 














9c = 7r 














0c = arccos (uja/"/'^) 


/ , ,2 \ 1/2 

eos0 


^(l-^)^/%m0 







C. Expectation Values of Field and Matter Observables 

As mentioned in Section A, the trial state is a tensor product of the coherent states 
of the one mode electromagnetic field times the collective atomic state. The expectation 
values of the quadratures of the field, g, p, and their fluctuations, can therefore be 
written in terms of real variables q and p. In a similar form, we can determine the 
expectation values of the photon number operator and its corresponding fluctuations. 
For the matter observables we have a similar situation, and the expectation values of the 
occupation operator Jz and atom dipole moments Jx and Jy can be expressed in terms 
of the stereographic projection variables {9, 0) indicating a point on the Bloch sphere. 

Then the expectation values of the quadratures of the electromagnetic field at the 
mimimum critical points, given in Eqs. lpfTOl) . are shown in Table ([Tll), where we have 
divided the expectation values by the appropiate normalization quantities. The fluc- 
tuations Ap and Ag are those associated to the coherent state of the electromagnetic 
field so they have their minimum values allowed by the Heisenberg uncertainty rela- 
tions, i.e., Ap = Ag = \/\/2. The expectation value of the number of photons and 
its corresponding fluctuation at the minimum critical points are also indicated in Table 
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TABLE III: Fluctuations of the operators J^, Jy and Jz divided by the square of the number 
of atoms in the system, at the minima. 



Minima 


(AJ,)7iV2 


(AJ,)ViV2 


(AJ,)ViV2 


0c = 


1 

AN 


1 

m 





ec = ^ 


1 


1 

AN 





Be = arccos 






47V 



The expectation values for J^, Jy, and Jz per particle are given by the spherical 
coordinates of a point on a sphere of radius 1/2. Therefore at the North and South 
Poles we get {Jx) = {Jy) = 0, and {Jz) /N = ±1/2, respectively. For the Parallel region 
we have 

{Jx) I uj\ {Jy) 1 ^A-, {Jz) i,^. 

= -^1 T-COS0, — ^ = -4/1 ^smffl, = . 13 

AT 2y 7^ ^' 2y 74 ^' AT 272 ^ ' 

The corresponding fluctuations per particle are indicated in Table fillip , their values tend 
to zero as the number of atoms in the system increases. 



Squeezing and Entanglement Entropy Se 



It is well known that coherent states for the electromagnetic field minimize the un- 
certainty relations. Radiation is said to be squeezed if the uncertainty of one quadrature 
is less than the standard quantum limit js, 21|. For spin or angular momentum systems. 



ever since the contribution by Kitagawa and Ueda [22| it has been recognized that the 
phenomena of squeezing is present when the fluctuations of the spin are correlated. It is 
then straightforward to prove that the SU (2) coherent states or collective atomic states 
are not squeezed and they have a squeezing parameter equal to unity. As a consequence. 
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the squeezing parameter 



3 

where J±_ denotes a component of J orthogonal to (J), gives us information of how good 
the trial state is approximating the exact solution of the eigenvalue problem of the TC 
Hamiltonian, i.e., the approximation is as good as how close the squeezing parameter is 
to one. 

Taking the partial trace for our trial state in the electromagnetic or in the atomic 
parts, there is no entanglement entropy. However, there is another way to separate the 
system and take the partial trace, by considering the Jordan-Schwinger map 

Jz = \{hlh~h\hY ^+ = 4^1, J~ = h\h, (14) 
where the bosonic operators 6j, fej, i, j = 1, 2, satisfy the commutator algebra of two 
independent harmonic oscillators. Consider taking the partial trace over the electromag- 
netic and the bosonic operators associated to the excited two level system; the result 
will be the same as if we take the trace over the radiation field and over the ground state 



level of the density matrix of the coherent atomic state [24 1. 

The coherent atomic state in terms of states of a two dimensional harmonic oscillator 
can be written as 

„t^o V^^; (1 + icH^ 

where the state |2j — n2^n2) = \ni,n2), with rii denoting the number of bi bosons 
(number of atoms in the lowest state) and n2 the number of 62 bosons (number of atoms 
in excited level) , with the constraint that the total number of quanta rii + n2 must equal 
the total number of atoms = 2j in the system. 

The partial trace over the first oscillator of the density matrix associated to this 
atomic coherent state is 

^ /1\T\ 

'n)(n| . (15) 



.-T.,.)^E(:)- 



(i+icn 



14 



TABLE IV: Occupation probability and entanglement entropy of the ground state of the TC 
model by means of the tensorial product of coherent states as trial state. 



Minima 


Pn 


Se 


0c = 














0c = arccos [uJa/"/'^) 




- En=0 Pn HPn) 



Therefore, the reduced density matrix is diagonaL Using Eq.([7]) the entanglement en- 
tropy is 

N 

SE{0) = -Y.Pni(^)HPn{e)) , (16) 



n=0 

where 



Pnie) 



N\ /l-cos^\"/l + cos^^ ^-"^ 



n J \ 2 

is the probability of a binomial distribution, which also corresponds to the occupation 
probabilty of n particles in the excited state of a two-level system. Through this expres- 
sion it is immediate to determine the estimates on the occupation probability and the 
entanglement entropy of the ground state. The results are given in Table (II VI) . 



D. Trial state in control parameter space 

The expression for the trial state \a, () that minimizes the energy surface takes 
the following forms: At the North Pole, uja > 7^, the ground state is given by 

li^np) = |0) ® |j, -j) . (17) 

At the South Pole, ua < —7^, the trial state has the form 

li^sp) = |0) ® |j, j) . (18) 
15 



For the Parallel case, \u!a\ < 7^, the approximate expression to the ground state is given 

by 

+j +00 

\lppar) ^ Y Yl ^ra,A^) ® |j, m) , (19) 



m=—j i/=Q 

where we define the expansion coefficients 



One usually would perform this sum by selecting a value of v to make the sum over 
m, and then proceeding with the following value of i^, until we reach some type of 
convergence. In the TC model A = m + is a conserved quantity. By replacing m by 
A — I/, we can write 

+3 ^+3 +00 A+j 

IV'par) = Y I]^A-.,.|J^) ® |j, A - I/) + ^ ^ Aa-.,.|z^) » |j, A - Z/) , (21) 
A=-j i/=0 A=j+1 i/=A-j 

which, in compact form, is 

+00 A+j 

IV'par-) = 5Z J] ^A-.,. «) |j, A - i/) . (22) 

A=— J i/=max(0,A— 

III. QUANTUM ANALYSIS 

A. Ground State Energies and Constant of Motion 

To understand the physical meaning of the constant of motion A = a^a + J^, we apply 
the transformation 

W(?7) =exp(i?7A) . (23) 
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to the TC Hamiltonian. The effect of the transformation in the quadratures of the 
electromagnetic field is 

('], (24) 
\p I \ — sin cos rj I \ p I 

making a rotation by an angle t] counterclockwise in the plane of the electromagnetic 
quadratures. In the same way we can calculate the induced transformation for the matter 
observables and Jy 

U(r,)('']uHr,) J (''] . (25) 

\Jy J ysm?7 COS?7 J \Jy J 

which is a rotation by an angle rj clockwise along the 2;-axis. This transformation then 
leaves invariant each term in the TC Hamiltonian. Therefore, we propose that the 
quantum states for the TC Hamiltonian have the form 

\^k)= c('=)|z/)®|j, A-z/) , (26) 

i/=max[0,A— 

because only one value of A is present. This contrasts with the trial function that we 
are considering, which is a superposition of all A's. The index k is denoting the different 
eigenvectors of the Hamiltonian matrix, whose dimension depends on the value of A: for 
X < j one has matrices of dimension d = X + j + 1, while if A > j the matrices have 
dimensions d = 2j + 1. We are using the natural basis states to diagonalize the TC 
Hamiltonian because it has the constant of motion A; for the general case, when one is 
not considering the rotating wave approximation, a new technique has been proposed 
to obtain the exact result for systems with size two orders of magnitude higher than 



reported in the literature [231]. 



In this contribution we will mainly be concerned with the behaviour of the ground 
state, which will be denoted by \ipgs) and its expansion coefficient by c^,. 
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Finding analytic solutions for small values of A is straightforward, for an arbitrary 
number of particles A^, any value of the detuning parameter A, and an arbitrary coupling 

r 

i) For A = — y we have a 1 x 1 Hamiltonian matrix, the ground state energy is given 
by Eo/N = - A), and the state is |0) O |j, -j). 

ii) For A = — Y + 1 we have a 2 x 2 Hamiltonian matrix, the ground state energy is 



given hj Ei/N = ^{2 - N - A + NA - a/It^TA^), and the state is a linear 
combination involving and 1 photons: |0) ® \ j, —j + 1) and |1) ® \ j, —j). 

iii) For A = — y + 2 we have a 3 x 3 Hamiltonian matrix, the ground state energy is 
given by 



, NA-N-2A + 4 2^M_^ + A2 sin (1(20 + vr)) 

""^/"^ = m ■ 

where N > 2 and 



tan ^ z= 

3v^72A 

V / 

and the state is a linear combination involving 0, 1 and 2 photons: |0) ® \j, —j + 2), 
|j, + and |2)® |j, -j). 

The expressions get more and more complicated, but one can solve analytically for 
up to A = — Y + 4. The way in which our Hamiltonian is written allows us to take the 
thermodynamic limit ^ oo in the expressions for the energy; in all cases mentioned 
these reduce to Eq/N. When in resonance A = 0, the results for all these values of A 



coincide with those given by Buzek et al. 



In Fig. ([5]) we show the quantum phase transitions for a system of 6 atoms (j = 3), 
with a detuning parameter A = 0.2. The straight lines correspond to the energy of 



the ground state for different values of the constant of motion A, starting from —3 (0 
photons) to +1 (4 photons). In each case the eigenstate is a mixture of to A + j 
photons. 



-0.40 




-0.45 - 



-0.50 - 



-O.55I1 



0.8 0.9 1.0 l.I 1.2 1.3 

FIG. 5: Ground state energy as a function of coupling constant 7. Results are shown for values 
of A from -3 to +1, for 6 atoms and A = 0.2. The heavy hne represents the minimum energy 
for all 7, and each change in A reflects a quantum phase transition. 



By taking the trace with respect to the field (matter) states the reduced density 
matrix takes the form 



min{A+f ,Af} 



^matter 



Q 



n=0 



field 



— n, n) {N — n, n\ 



where c^+n_^ 



or c,y) is determined from the Hamiltonian diagonalization. Taking the 
trace over the first oscillator in -we have a diagonal matrix of the same size given 

by 



^A+f-ni 



■'ni,n2 



(27) 
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where ni = 0, ■ ■ • , min{A + y, A^}. These expressions will be used below to calculate the 
fidelity between the variational and the exact quantum states. As the reduced density 
matrix of the matter is diagonal, the matter-field entanglement entropy equals that 
between the atoms occupying the two hyperfine levels 

min{X+^,N} 
n=0 

in both cases. 

The expectation values for the photon, Ojr, and atomic, observables with respect 
to the quantum ground state can be simplified to the expressions 

A+i 

(^g^ldpl^g,) = \Cuf {'^\OfW) (28) 

f=max[0,A— j] 

A+i 

Oa \^gs) = Yl l^-l' \j, A - z/) . (29) 

i/=max[0,A— j] 

Given the form of the quantum ground state ( l26i) it is immediate that: 1) the expec- 
tation values of any power of the ladder operators J± vanish because A and the number 
of excited atoms cannot be changed when the number of photons remains constant. 
Therefore (J^) = {Jy) = 0, and (J^) = (J|) because we have only contributions from 
the diagonal part of the operators. 2) The expectation values of any power of the cre- 
ation and annihilation operators vanish for the same reason. In particular we obtain 
(g) = {p) = {{qp + pq)) = 0. For the expectation values of and p^, one has again 
contributions from the diagonal part and they are given by 

where u denotes the number of photons of the state. 3) The fiuctuations in the quadra- 
tures of the fiels are given by (Ap)^ = (g^) and (Ag)^ = {p"^). 

As in the TC Hamiltonian the expectation value of the angular momentum vector is 
only in direction z, and the fiuctuations of the dipole moment components of the atoms 
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satisfy (AJ^)^ = (AJ^)^; the squeezing spin coefficients in directions x and y must be 
equal. Therefore, the expression for ^, in the orthogonal directions to the z-axis, is 
defined by the expression 



where k = x or y. 

IV. SEMICLASSICAL VERSUS QUANTUM RESULTS 

A. Ground State Energies and Constant of Motion Expectation Values 

Fig-Q displays the energy per number of particles Eq/N, and the constant of motion 
Ao, as a function of the coupling interaction 7. The top figures correspond to a number 
of atoms = 6, while the bottom figures correspond to = 100; in both cases the 
detuning parameter is A = 0. While the classical energy slightly overestimates the 
exact quantum result outside the North Pole (as must be according to the Rayleigh-Ritz 
variational principle), second order corrections to the eigenvalues of the Hamiltonian 
matrix gives a much better agreement when the number of atoms is increased, as shown. 

We observe that both Eq and Aq depend strongly on the coupling strength 7. As 
this coupling is increased, the value of the (quantum) constant of motion at which the 
minimum energy is obtained remains constant until a threshold is reached, at which 
point it jumps to its next value. This same process repeats itself as 7 keeps increasing, 
Ao taking discrete values throughout (in contrast to its semiclassical counterpart, which 
is continuous), reflecting the quantum nature of the system. As increases, the width of 
the steps becomes smaller until the difference with the classical result seems negligible. 
However, the separation of A- values always equals 1 (cf. zoom in the figure). In the 
North Pole, the state is |z/) ® \j, m) = |0) (g) \j, —j), so only one value of Aq (and of Eq) 
is had all through this region — 1 < 7 < 1, and it naturally coincides with the classical 




(30) 
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FIG. 6: The energy per number of atoms, Eq/N (left), and the constant of motion, Aq (right), 
are shown as functions of the interaction parameter 7 for both, the semiclassical and quantum 
cases. The exact quantum energy is equal to or less than its classical counterpart. The 
semiclassical constant Aq has continuous values while the exact quantum result is discrete. We 
use = 6 (top), and N = 100 (bottom). In both cases A = 0. For N = 100 the graphs totally 
superimpose; however, the zoom in the last graph makes the discreteness in Aq evident. 

value. The same (dual) result is obtained in the South Pole. 

This discrete behaviour in Aq is inherited in all the quantum observables of interest, 
as we shall see below. It is important to note the horizontal shape of the steps in Aq 
shown in Fig.([HD. This is a consequence of the system being in resonance: A = 0. 
When the detuning parameter is taken away from resonance, the Aq steps acquire a 
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slope given precisely by the difference in energy between the absorbed photons and the 
atom's energy level separation. 

Given the present possibility of confining a small number of atoms and performing 
non-demolition experiments on them, and the direct importance of these on cavity QED, 
quantum information theory, encryption, teleportation, and quantum optics in general, 
in all the results that follow, we take = 6 and A = 0.2 in order to illustrate the 
behaviour of a small system slightly away from resonance. The effect of the detuning 
parameter on the ground state energy, for A = —0.2, 0, 0.2, is shown in Fig. ([7]). One 
notices that the width of the North Pole depends on the value of A, as expected. 
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FIG. 7: Energy per number of atoms, Eq/N as a function of the parameter 7. We use j = 10. 
The lowest energy curve corresponds to A = —0.2, the intermediate energy curve to A = 0, 
and the curve with larger values to A = 0.2. 

B. Matter Observables 

Fig® presents the expectation value of {Jz)/N for = 6 atoms and A = 0.2. The 
observable takes discrete values for the exact quantum solution while it is continuous for 



23 



the semiclassical trial state. The inherited discreteness from Aq, associated to quantum 
phase transitions as the structure of the ground state changes with 7, is evident. The 
semiclassical approximation follows the same trend. 

For 7^ < uja all the atoms are concentrated in their lowest hyperfine level. It is natural 
to find that {Jz)/N takes on the value —j/N, which in the plotted case corresponds to 
—0.5, and the constant of motion takes its minimum value A = —3. For values of 7 
outside this region there are other occupancies. When 7 is very large the occupancy 
tends to zero, i.e., the two hyperfine levels are equally occupied. The dispersion (AJ^)^ 
has values ten times smaller. 
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FIG. 8: The expectation value {Jz)/N (left) for the semiclassical (continuous) and quantum 
(discrete) models, together with their corresponding fluctuations (right), are displayed as func- 
tions of 7, for = 6 and A = 0.2. 

C. Field Observables 

The expectation value {n)/N of the photon number operator per number of atoms 
with respect to the ground state is shown in Fig. ([9]), for the exact and variational 
quantum results, for a system of = 6 atoms with a detuning A = 0.2. The fluctuation 
in the number of photons {An)'^/N'^ is also shown. Note that in the North Pole, {u) ® 
\j, m) = |0) ® |j, —j), the number of photons is zero. Thus, the fluctuation must be 
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zero and completely coincide with its semi classical counterpart. Outside this region, the 
proposed trial function as a coherent state for the description of the electromagnetic field, 
which includes an infinite number of photons in its composition, greatly overestimates 
the quantum result: for 7^2 there is approximately one photon per atom with very 
good agreement between the two procedures, while their corresponding fluctuations have 
at least a difference of a factor of six. 
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FIG. 9: The expectation value of {n)/N (left) and its corresponding fluctuation (An)^/Af^ 
(right), are shown as functions of 7. We use = 6, and A = 0.2. 



D. Squeezing, Entanglement and Fidelity 

We have found that the ground state of the TC model does not present the phe- 
nomenon of squeezing in the matter components. However we propose to use the be- 
haviour of the squeezing spin coefficient ^ as an additional signature of the goodness 
of our trial state in reproducing the exact quantum solution. We can see in Fig. lfTOl) 
that the squeezing coefficient exceeds the value of one for large values of the coupling 
parameter. This means that the trial state must have a very small overlap with the exact 
quantum state, and it is due to the multiple values of the constant of motion contained 
in the composition of our trial state. 

The entanglement is usually regarded as a purely quantum correlation which plays 
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an important role in quantum phase transitions 10|, lUj • For the variational state the 
entanglement entropy between the matter and field is zero because it corresponds to a 
pure state (matter coherent state). It is then interesting to do a comparison with the 
entanglement between atoms in the semi-classical approach. Fig. lITOl) shows the entan- 
glement entropy associated to the variational and exact quantum cases, for a reduced 
density matrix that eliminates the electromagnetic degrees of freedom and discriminates 
the matter part into the atoms occupying the lowest hyperfine level from the rest, and 
viceversa. Once again, the semiclassical result is a very good approximation to the exact 
one. The plots shown are for = 6 and A = 0.2. It is clear that for any A > — j (i.e.. 



outside I7I < y/1 — A) there is entanglement between the atoms. In the thermodynamic 
limit N 00, the behavior is the same with an even more abrupt slope at the separa- 
trix. Besides, one can notice directly from this plot, that the entanglement entropy with 
7 ^ 2 is S'^; ~ 1.6 nats; this compares favorably with a straightforward estimate of the 
maximum entanglement entropy available for a quantum system of j = 3, which would 
yield Se = 1-94 nats. 
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FIG. 10: The entanglement entropy (left) and the squeezing spin coefficient (right) are plotted 
as functions of 7, for = 6, and A = 0.2. 



An interesting measure, coming from quantum information theory, is the fidelity 
F{pii P2) 25|, where pi denotes the density matrix of a system. This is a measure of 
how close two probability distributions are to each other, even for mixed states. In the 
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case of two pure states, the measure is equivalent to the overlap. In our case, we want to 
compare the reduced probability distributions calculated from the exact model and from 
the variational model. Since they are both diagonal (cf. Eqs.ffTSl [271)). the expression 
for the fidelity is simple: 

min{A+f ,Af} 

Fig^',go')= Yl 7|cA+f-nl>n(^c) • (31) 
n=0 

Fig. (fTT|) shows the fidelity as a function of 7 both in resonance A = 0, and away from 
resonance A = 0.2, for the probability distributions mentioned above. At the North 
Pole (zero photons and all atoms in the ground state, A = — j) the fidelity represents 
the overlap between the two states, and equals 1. Each jump represents a value of A 
increased by 1. The first and largest transition is found at the separatrix, while the other 
discontinuities occur at the quantum phase transitions. Note that F diminishes at each 
of these for small values of A (small number of photons). In all these cases we have mixed 
states. After X = j (number of photons equal to the number of atoms) the dimension of 
the Hamiltonian matrix remains constant and F increases at each transition (except for 
the anomalous case = 1). It is important to stress the fact that the fidelity is high 
and the probability distribution of the coherent state approximates very well the exact 
solution, especially for a small number of atoms and a large number of photons. 

E. Occupation Probability Distributions 

The composition of the ground state of the TC model can be analyzed both, in terms 
of the distribution of the atoms into the two hyperfine levels, and of the photon number 
content, be it for the exact energy ground state or for its variational approximation. 
Fig. (1121) shows the composition of the trial and exact states for A^ = 6, A = 0.2. The 
four plots show the change in the occupation probability as one crosses one arm of the 
separatrix: in this case, we move along Path II in Fig.([T]) and cross the separatrix at 
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FIG. 11: Fidelity as a function of 7, for j = 3. Light curves correspond to resonance A = 0, 
and dark curves to A = 0.2. 

7 = — \/0.8 from the Parallels region into the North Pole (recall that A = 0.2). The 
marked change due to the phase transition is obvious. While the composition of the 
ground state in both cases are very similar, one has to notice that the variational state 
is constituted by a distribution of A's in contrast to the exact quantum result which 
is described by only one A-value. (This good approximation is a consequence of the 
fact that the distribution of A's is a gaussian centered at the quantum eigenvalue, as 
will be seen below in Fig.llMl) and its corresponding discussion.) This quantum A-value 
depends of the coupling parameter 7 considered in the TC Hamiltonian as can be seen 
in Fig. ([6]). It is important to stress that, for values of 7 far away from the separatrix 
(7 = —1.5 in the figure), the semiclassical occupation probability estimates very well the 
quantum result. The estimate is even better as N increases, until they become almost 
indistinguishable, as shown in Fig. ([12]) calculated for = 100. However, as we get 
close to the separatrix (7 = — -s/cTS for = 6; 7 = — 1 for = 100) the semiclassical 
approximation gets poorer. When we cross it, both distributions become a spike at 
X = m = —j , and they remain so until the separatrix is crossed again. The distributions 
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are exactly the 



same changing 7 for —7 (cf. Fig. fll3p ). 
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FIG. 12: The composition of the trial and exact states of the TC model are displayed for = 6, 
A = 0.2, as we cross the separatrix along Path II for values of 7 = —1.5, —1.0, —0.9, —0.8. 
The corresponding quantum eigenvalues of the constant of motion are A = 2, —1, —2, —3. m 
runs from —j to j (—3 to 3 in this case) along the horizontal axis. The narrow bars correspond 
to quantum values. 

To calculate the distribution of the A values for the trial state one can use the coeffi- 
cients fl20p . Fig. (1141) shows the obtained distribution for 7 = —1.5, = 6, and A = 0.2. 
Notice that there are approximately six different values of A with significative values 
in the probability distribution. The maximum is centered between A = 1 and A = 2; 
adjusting a gaussian distribution we get a mean /i = 1.87 and a standard deviation 
a = 2.06. For reference, the exact diagonalization gives Aq = 2 for this value of 7. While 
only values of A up to 10 were considered in this chart, these amount to 99.979% of the 
total contribution. This is the main reason why the tensorial product of coherent states 
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FIG. 13: The composition of the trial and exact states of the TC model are displayed 
for N = 100, A = 0, as we cross the separatrix along Path II for values of 7 = 
—2, —1.5, —1.1, —1.01, —0.9, 1.01. The corresponding eigenvalues of the constant of motion 
are A = 81, 23, —31, —48, —50, —48. m runs along the horizontal axis. The narrow bars 
correspond to quantum values. 



constitutes a good approximation to the quantum ground state. 

In order to explain the good approximation of the trial state in all the studied ob- 
servables except for the fluctuations in the number of photons, one can study the way in 
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FIG. 14: A distribution of the variational state, for = 6 atoms, A = 0.2, and 7 = —1.5. It is 
approximately a gaussian distribution with mean /i = 1.87 and standard deviation a = 2.06. 
The exact diagonalization gives Aq = 2 for this value of 7. 

which the number of photons is distributed in the ground state for each quantum eigen- 
value A. From our analysis on fidelity (cf. Fig. lfTTl) ). we expect the distributions for the 
variational state and the quantum state to agree inside the North Pole and away from 
the separatrix, but to show a strong disagreement at the separatrix itself. This is exactly 
what Fig. (fT5|) shows, where, starting at the minimum value A = —10, corresponding to 
7 = a/1 — a — 0.01 (just inside the North Pole), we increase A in steps of 1, crossing the 
separatrix and reaching a value away from it at A = 6, corresponding to 7 = 1.5. We 
see that, even though the photon number distribution of the trial state is wider than 
the quantum case, the mean is the same. The top of Fig. fll6l) . for = 10, A = 0.2 and 
7 = 5 (corresponding to A = 124), shows this more clearly. Considering the contribution 
to the photon number probability distribution given by only the A value of the exact 
quantum solution, we obtain the narrow, darker bars in Figs.(fT5| [T6l) . The amplitudes 
are necessary smaller since we have eliminated all contributions from other A values: 
for instance, with = 10, A = 0.2 and 7 = 5 the contribution of the values of v to 
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the state is of only 3.5%. However, a proper renormalization renders this distribution 
practically equal to that of the corresponding quantum photon number distribution, as 
illustrated in Fig. f|T6|) (bottom). 
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FIG. 15: {color online) Photon number distribution for different eigenvalues of the constant of 
motion A, for the exact quantum ground state with = 20 particles and A = 0.2, represented 
by the medium-width (blue) bars. Starting at the minimum value A = —10, corresponding 
to 7 = Vl — ^ — 0.01 (just inside the North Pole), we increase A in steps of 1, crossing the 
separatrix and reaching a value away from it at A = 6, corresponding to 7 = 1.5. In wider bars 
(magenta) the corresponding distribution for the variational case, containing all values of A, 
is shown. If we consider the trial state with the contribution of only the A value given by the 
exact solution, the narrow (darker, green) bars result. 
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FIG. 16: Above: Same as Fig. ijlSp for A = 124, corresponding to 7 = 5, very far away from 
the separatrix. Below: Comparison of the renormalized trial state distribution (narrow bars) 
with the quantum distribution, for the same values. 
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V. CONCLUSIONS 



A comprehensive study of the ground state quantum phase transitions in the Tavis- 
Cummings model, for any number of atoms, has been presented as a function of the 
strength of the interaction 7 and the atom hyperfine two- level separation energy cua- 
The minima of the expectation value of the Hamiltonian were calculated for a coherent 
trial state, the separatrix that divides the control parameter space was found, and the 
structure of the ground state of the system in all these regions was studied. 

First and second order phase transitions were determined as the separatrix is crossed, 
and their influence on the behaviour of observables of interest for the matter and the field 
were explored. These include the population and dipole moments of the two level atoms, 
the expectation value of the number of photons and its fluctuations, and the squeezing 
and entanglement properties of the ground state. Additional quantum phase transitions 
of the ground state of the model, due to the discreteness of the value of the constant 
of motion Aq, disappear when the number of atoms grows to infinity. However, the 
expectation values of the field and matter observables, together with the entanglement 
entropy, are manifestly discrete even for a very large number of atoms. 

A discrete behaviour in the expectation value of the constant of motion A at which 
the minimum energy is obtained, reflecting the quantum nature of the system, was 
shown: as the matter-field coupling is increased, the value of Aq remains constant until 
a threshold is reached, at which point it jumps to its next value. The discrete behaviour 
is inherited in all other quantum observables of interest. Since the value of observables 
remain constant over small regions of the interaction parameter, this feature can be used 
for control and manipulation of quantum systems of small number of atoms, where non- 
demolition experiments could be carried out. It was also found that when the detuning 
parameter is taken away from resonance, the Aq steps (and those of other observables) 
acquire a slope given precisely by the difference in energy between the absorbed photons 
and the atom's energy level separation. The effect of the detuning is of an even greater 
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importance, as it modifies expectation values of interest sucli as tlie number of photons, 
the entanglement entropy, and the dipole moments. Our calculations show that these 
increase as A decreases, and viceversa. Thus, one could fine-tune the entanglement 
entropy Se between atoms, for instance, through a fine-tuning of A. 

While their semiclassical counterparts show, in contrast, a continuous behaviour, we 
have shown that they are an excellent approximation to the exact (quantum) values in 
all the calculated observables, with the exception of the photon number fluctuation. The 
reason for this is found in the fact that the proposed trial function, as a coherent state 
for the description of the electromagnetic field, includes an infinite number of photons 
in its composition, greatly overestimating the quantum result. This results in the trial 
state having a very small overlap with the exact quantum state, which may also be seen 
in the squeezing coefficient, as it exceeds the value of one for large values of the coupling 
parameter (cf. Fig. fllOp ). In reality, the variational state has a gaussian distribution 
of X's around the quantum value of the constant of motion. In the case of = 6 the 
standard deviation of this distribution is cr = 2, implying that significant contributions 
will be had only for a small number of occupation states, or equivalently, from a small 
number of Fock states. A better trial state may be constructed by truncating the 
proposed state to only one value of A, which can be chosen to be in agreement with that 
of the quantum solution, and by renormalizing the state. When this is done, we have 
found an excellent agreement between this renormalized trial state and the quantum 
result (cf. Fig. (fT6l) ). The calculation of the expected values of other observables, using 
the renormalized state, is undergoing. 

The composition of the ground state was shown in terms of the distribution of the 
atoms into the two hyperfine levels, both for the exact energy ground state and for 
its variational approximation. We showed that away from the separatrix the semiclas- 
sical occupation probability estimates very well the quantum result. The estimate is 
even better as the number of atoms increases, until they become indistinguishable. 
However, as we get close to the separatrix the semiclassical approximation gets poorer. 
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The fidelity of the reduced probabihty distributions of the atoms populating the ex- 
cited hyperfine level is calculated for the semiclassical and exact ground states, which 
gives information about the overlap in the region where a pure state dominates, i.e., 
I^aI ^ 7^- Outside this region we obtain information about the similarity between two 
mixed density matrices. This fidelity can also be fine-tuned through A. 
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